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hancement of low-m,t kaons (K^ and K ) in Si+Pb collisions at AGS can be 
explained if a density isomer is introduced in the nuclear equation-of-state. 
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I. INTRODUCTION 

An interesting experimental observation in relativistic heavy-ion collisions at AGS en- 
ergies is the enhancement of low-rrit kaons (both K~^ and K^) in Si+Pb collisions at 14.6 
(GeV/c)/nucleon |l|. The spectra of these extremely 'cold' (in terms of their transverse 
mass rrit) kaons can be characterized by an inverse slope parameter (temperature) as low as 
15 MeV, which is about one order-of-magnitude smaller than the temperature of 'normal' 
kaons measured in earlier AGS experiments 0]. These low rrit kaons cannot be obtained 
in conventional calculations such as the relativistic quantum molecular dynamics |^. At 
lower incident energies, a similar enhancement of low-pt pions has already been observed at 
Bevalac [Q and SIS energies [^, and one of the plausible explanations for this enhancement, 
as put forward by Xiong et al. , is based on the attractive pion potential due to the A- hole 
polarization. 

Very recently, Koch has tried to explain the enhancement of low-r/it kaons at AGS 
energies by a similar mechanism as in Ref. 0, i.e., the attractive kaon potential. When a 
kaon 'climbs' out of a potential well, it has to lose its kinetic energy and thus gets cooled 
down. It was found in Ref. [0], however, that in order for the attractive potential to be 
effective in cooling down the kaons, the baryon fireball, which is the source for the kaon 
potential, has to expand slower than kaons, so that the latter will have the chance to climb 
out the potential well before it disappears. 



In a relativistic transport model fg,^ based on the a-uj model |I0| for the nucleon- 
nucleon interaction, as the one used in this work and in Ref. [0, the expansion of the fireball 
is controlled by the delicate balance (or rather unbalance) between the scalar (a) attraction 
and vector {uj) and thermal repulsion. In usual relativistic transport model, the scalar and 
vector coupling constants, which are determined by fitting the properties (saturation density, 
binding energy, effective mass and incompressibility) of normal nuclear matter, are assumed 
to be density and temperature independent. For a fireball at about six times normal nuclear 
matter density and a temperature of about 170 MeV, as expected for the AGS energies 



PJTTH, the scalar attraction is overwhelmed by the vector and thermal repulsion, as the 
former tends to saturate at high density due to relativistic effects. The fireball would thus 
expand rather fast and kaons cannot be cooled down effectively. 

To slow down the expansion of the fireball, Koch ||^ has suggested that the system 
initially goes through the chiral restoration transition, which leads to a softening of the 
equation-of-state and thus slows down the expansion of the fireball. This has been achieved 
in Ref. by reducing the vector coupling constant to about half of the value used in the 
original Walecka model |T^. Using the coupling constants given in Ref. [^, however, the 



nuclear matter would have saturated around 0.6 fm~^ with a binding energy of about 250 
MeV, which is clearly unrealistic. 

In this paper, we propose that the slowing down of the fireball expansion is also possible 
if there is a density isomer at high densities and if the fireball is initially at this abnormal 
state. Lee and Wick [|1^] observed that the chiral sigma model can lead to an abnormal state 
at high densities. The binding energy of this state can be very large, leading to a secondary 
minimum in the nuclear matter equation-of-state. The chiral sigma model, however, cannot 
describe correctly the saturation properties of the nuclear matter. In a relativistic mean- 
field theory that includes the A degree of freedom, Boguta |TB[ has shown that it is possible 



to simultaneously describe the nuclear matter properties at saturation density and predict 
the existence of a density isomer at high density. This has been realised in Ref. |]13| by 
assuming that the scalar interaction of the A is stronger than that of the nucleon. In the 
present exploratory study, we will simulate the effect of the density isomer by introducing a 
density-dependent vector coupling constant in the non-linear a-u model, which gives correct 
saturation properties of the nuclear matter at po = 0.16 fm~^ and leads to a secondary 
minimum around (5-6)po- More consistently, one would like to combine the chiral sigma 
model with the Walecka model, so that both the nuclear matter saturation and the density 
isomer come out naturally from the underlying Lagrangian. 



II. NUCLEAR EQUATION-OF-STATE 

Our study is based on the relativistic transport model for the expansional stage of the 
fireball formed in a Si+Pb collision at 14.6 (GeV/c)/nucleon |^^15|. In Ref. |16|, we have 



proposed two sets of parameters for the non-linear a-uo model. In this work, we use slightly 
different parameters which are determined by requiring a saturation density of 0.16 fm~^, a 
binding energy of 15.96 MeV, a nucleon effective mass m* /m=0.77, and an incompressibility 
of 380 MeV. Using the notation of Ref. [|I^, these parameters have the following values 

Cv = 10.5, Cs = 13.61, B = -0.003131, C = 0.02908. 

The nuclear equation-of-state corresponding to this parameter set is shown in Fig. 1 by the 
dashed curve. To simulate the density isomer at high densities, we assume that the vector 
coupling constant is density-dependent and is determined by nuclear matter properties at 
saturation density and the assumption that around 5.5po a secondary minimum appears. 
Explicitly, we have 

Cv = 10.5, p < po, 

= 10.59 - 0m23(—Y, po<p< 5.5po, 

= 7.8, p>5.5po. (1) 

The nuclear equation-of-state with this density-dependent vector coupling constant is plotted 
in Fig. 1 by the solid curve. A secondary minimum is seen around 5.5po with a binding 
energy of about 2 MeV. 

We note that the density dependence of the vector coupling constant leads to a so-called 
rearrangement potential for the nucleon. In the present case, this potential is attractive, so 
it tends to further slow down the expansion of the fireball. 

The fireball, which contains nucleons (A^), deltas (A), hyperons {Y), pions (tt), rhos (p), 
as well as kaons (K^) and antikaons (K^), is assumed to have an initial density of 6po and 
temperature of 170 MeV. The initial density is higher than that used in the fireball model of 



Ref. |T^ but is consistent with that predicted by the cascade model |T^ and the relativistic 



quantum molecular dynamics . The motions of these particles are then described by the 
relativistic transport model. The treatment of two-body scattering is the same as in Refs. 



I5| , p!6[] . For kaon-baryon and antikaon-baryon elastic scattering we use the empirical cross 
sections [^ . The kaon-pion and anitkaon-pion scattering are dominated by K* resonances, 
and we use the Breit-Wigner cross section of Ref. |jl9[|. For an antikaon, it can also be 
annihilated by a nucleon into a hyperon and a pion {K~N -^ Ytt). The annihilation cross 
section is appreciable and is expected to be modified strongly in a medium due to changing 
in-medium hadron masses. The in-medium cross section for antikaon annihilation will be 
discussed later. 

III. KAON AND ANTIKAON POTENTIALS 

In a nuclear medium, a kaon feels an attractive scalar potential due to explicit chiral sym- 
metry breaking and a repulsive vector potential due to tu-exchange. The vector potential for 
an antikaon is attractive under the G-parity transformation. In mean-field approximation, 
the kaon and antikaon dispersion relations can be written as 

*2 2 , 1 2 '^KN I 3 CU^' 

Ik ^ Ik 

2 , 1 2 ^KN , 2 w 5f^ \2 

= m^ + k^ - -T^ps ± i:^k[ ) Pb, (2) 

Ik "J "'■1^ 

where ps and ps are the baryon and scalar densities, respectively. The second line follows 

from the SU(3) relation g^^ = 3gp and the KFSR relation nip = 2\/2fgp (assumption of 

/i^ = A = / is implied) 0. 

In Eq. (2) the plus and minus signs are for kaon and antikaon, respectively, T^kn is 

the KN sigma term, and fx is the kaon decay constant. In this work, we use T^kn = 350 

MeV and fK= 93 MeV, as used in our study of kaon and antikaon production in heavy- ion 



collisions at SIS energies [^^ . Following Ref. p3|, we define the kaon (antikaon) potential 
as 



UK,K{P^k)=UJ*{p,k)-COo{k). (3) 

The results for the kaon (antikaon) potential with zero momentum in a nuclear medium 
is shown in Fig. 1. One sees that the kaon potential is weakly repulsive at low densities, 
qualitatively consistent with that from the impulse approximation based on the kaon-nucleon 
(KN) scattering length. At high densities (above about 3po) the kaon potential becomes 
attractive as a result of the reduced repulsive vector interaction. The antikaon potential is 
attractive at all densities. 

The propagation of kaons and antikaons in their mean-field potentials is given by the 
following equations of motion 

* = 1^' li= -^^t'-^*(''-k). (4) 

r 2 1 1/2 

where E* = [m\ + k^ - {T^kn/PkIPs + {\^pBf\ , and Uk,k is given by Eq. (3). When 

a kaon (antikaon) propagates from the high density to the low density (especially when it 

crosses the surface of the fireball), its momentum is reduced by the attractive potential. 

This is the mechanism by which the attractive potential cools down particles. 

IV. ANTIKAON ANNIHILATION IN MEDIUM 

A main difference between a kaon and an antikaon in a dense matter is that the latter 
can be annihilated by a nucleon into a hyperon and a pion, i.e., K~N -^ Yn. This cross 
section at low energies has been analysed by the K- matrix method of Martin et al. [E^, and 
the results can be fitted by the following parameterization 

"""""■^"^'^ 1.0 + 76.7(v/i -^) + 161.1(v/i - ^f ' ^""^ 



where ^/s is the invariant energy of the K'N system and ^/sQ = rriN + tuk- In nuclear 
medium, we assume that the form of the annihilation cross section does not change, but 
replace a/s and ^/so by v^ and a/s^ = m^ + ^jv, respectively. Another important medium 
effect is that the reaction K~N -^ Yn may be forbidden in the medium. To see this, we take 
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the scalar potential of an antikaon to be approximately 1/3 of the nucleon scalar potential Us 
and that of hyperon to be {2/3)Us- In the reaction K~N -^ Yn, the total scalar potential 
is then (4/3) f/5 in the initial state and is more negative than that in the final state, which 
is (2/3) f/^. For antikaons with low momenta, which we are interested in, there is thus a 
possibility that the invariant energy ^/i* of the antikaon and nucleon is below my + m-jr, so 
the antikaon cannot be annihilated. 

V. RESULTS 

Because of the density isomer around 5.5po; the fireball expands rather slowly. In Fig. 3, 
we show the baryon density profile of the fireball at different times. In the first 4 fm/c, the 
central density decreases by only about 10%. From t=4 fm/c to 8 fm/c, the central density 
decreases from 5.5po to about 2.0po- Up to 10 fm/c, a well-defined surface is seen, and when 
kaons move across this surface their momenta are changed. 



To obtain good statistics, we use the perturbative test particle method of Ref. [^ for 
kaons and antikaons. Furthermore, we average the transverse mass spectra of kaons and 
antikaons over the entire rapidity range as in Ref. . Using a bin size of 5 MeV, we have 
more than one thousand events within each transverse mass bin. In Fig. 4, we show the 
transverse mass spectra of kaons at different times. The upper-left panel gives the initial kaon 
spectrum which are in an equilibrium fireball at temperature T = 170 MeV. The spectrum 
can be fitted by exp{ — mt/T') with an apparent temperature T' ~ 145MeV, which is smaller 
than the source temperature T. This difference is mainly due to the rapidity effect, as the 
apparent temperature T' is related to the source temperature T by T = T'cosh (y), where 
y is the kaon rapidity in the center-of-mass frame of the fireball. 

The upper-right panel gives the kaon spectrum after 4 fm/c. During this time interval, 
kaons with high momenta have mostly escaped the potential well. We already see an en- 
hancement of low-mt kaons and a weak low temperature component. The steep rise of the 
extremely low-mj kaons happens mainly during the time interval from t=4 fm/c to 8 fm/c. 



as shown in the lower-left panel of the figure. During this time interval, kaons that are 
initially trapped inside the potential well gradually climb out the well as the system slowly 
expands. A well-defined low temperature component is visible at this time. Afterwards 
there is little change in the kaon spectrum, as shown in the low-right panel for t=12 fm/c, 
as the density of the system is already below 2po and the kaon potential is very weak. 

In Fig. 5, we compare our results with preliminary experimental data of the E814 
collaboration |]l|. The experimental data are in the forward rapidity (2.2 > y < 2.4), while 
our results are averaged over the entire rapidity region as mentioned early. Since we do 
not have a complete dynamical model for the entire process of heavy-ion collisions at AGS 
energies, we scale our final spectrum such that we fit the data point at mt—niK ~ 7 MeV. The 
initial spectrum, given by the dotted histogram, is similarly scaled, and can be characterized 
by one apparent temperature of about 145 MeV. The final kaon spectrum, obtained after 
including both propagation in the mean-field potential and rescattering with baryons and 
pions, is given by the solid histogram and shows clearly a two-component structure. The 
component corresponding to high transverse masses {mt — mx > 0.02 GeV) has an apparent 
temperature of about 140 MeV, and is similar to the initial kaon spectrum. The steep rise 
of low-nit kaons (rrit — rriK < 0.02 GeV) can be fitted by an exponential with an apparent 
temperature of about 40 MeV. Although this temperature is still about a factor of two larger 
than the temperature extracted from the experimental data, it is considerably smaller than 
that of kaons with high transverse masses. 

Finally, we compare in Fig. 6 our results for antikaon spectra with preliminary data 
of Ref. H]. Again, we scale our final antikaon spectrum such that we fit the data point 
at rrit — tuk ~ 7 MeV. Because of annihilation, the final antikaon multiplicity is smaller 
than the initial one. For comparison, we have reduced the initial spectrum by the ratio of 
the final antikaon multiplicity to the initial one. The initial spectrum (dotted histogram) 
can be fitted by a single exponential with an apparent temperature of about 145 MeV. The 
final antikaon spectrum again shows a two-component structure. The high-rrit component is 
similar to the initial spectrum, while the \ow-mt component has a much lower temperature 
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of about 30 MeV, which is again about a factor of two larger than that extracted from the 
experimental data. The K~ temperature is lower than that of K^ as its potential remains 
attractive throughout the whole expansion of the fireball while the kaon potential becomes 
repulsive once the density is below about 3po- We note that the reduced annihilation of 
low-energy K^ in dense matter as discussed in section IV is crucial in obtaining the low-m^ 
antikaons. Finally, we mention that in our calculation the proton spectrum does not show 
a two-component structure, which is consistent with the experimental data. 

VI. CONCLUSIONS 

In summary, we have studied the enhancement of low-rrii kaons and antikaons in heavy- 
ion collisions at AGS energies in the relativistic transport model. By introducing a density 
isomer at high densities, the expansion of the fireball is slowed down so that kaons and 
antikaons can be effectively cooled down by their attractive mean-field potentials. We see 
clearly a 'cold' low-m^ component in both kaon and antikaon transverse mass spectra, as 
observed in AGS experiments by the E814 collaboration. To carry out more detail compar- 
isons with the experimental measurements, we need to introduce rapidity cuts in analyzing 
the kaon and antikaon spectra. Also, it is necessary to extend the relativistic transport 
model to include more degrees of freedom so that it will be suitable for describing the entire 
collision dynamics at AGS energies. 
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Figure captions: 

Fig. 1: Equation-of-state of nuclear matter. The dashed curve is based on the non- 
hnear a-u model parameters of Ref . ||T^ . The solid curve is based on the density-dependent 
vector coupling constant of Eq. (1). 

Fig. 2: Kaon and antikaon potentials (with zero momentum) in nuclear matter as a 
function of baryon density. 

Fig. 3: Baryon density profile of the fireball at different times. 

Fig. 4: Kaon transverse mass spectra at different times. 

Fig. 5: Initial and final kaon transverse spectra. The experimental data for 2.2 > y < 
2.4 are from Ref. [g. 

Fig. 6: Same as Fig. 4 for antikaons. 
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